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We study ground- and excited-state properties of all sd-shell nuclei with neutron and proton 
numbers 8 ^ N, Z ^ 20, based on a set of low-resolution two- and three-nucleon interactions that 
predict realistic saturation properties of nuclear matter. We focus on estimating the theoretical 
uncertainties due to variation of the resolution scale, the low-energy couplings, as well as from 
the many-body method. The experimental two-neutron and two-proton separation energies are 
reasonably well reproduced, with an uncertainty range of ~ 5 MeV. The first excited 2 + energies 
also show overall agreement, with a more narrow uncertainty range of ~ 500 keV. In most cases, 
this range is dominated by the uncertainties in the Hamiltonian. 

PACS numbers: 21.10.-k, 21.30.-x, 21.60.Cs, 27.30.+t 


Introduction. Recent advances in nuclear theory have 
established the importance of three-nucleon (3N) forces 
in understanding the structure of medium-mass nuclei, 
for the evolution to the neutron and proton driplines [1 
5] and the formation of shell structure [6-13]. Three- 
nucleon forces are also key for realistic saturation proper¬ 
ties of nuclear matter [14-18], which in turn are obtained 
from global analyses of all nuclei. To date, ab initio stud¬ 
ies of medium-mass nuclei have largely focused on closed- 
shell nuclei or isotopic chains, generally in the vicinity of 
semi-magic nuclei, and no comprehensive study exists to 
explore nuclear forces over a full range of the nuclear 
chart, such as the sd shell. 

An additional challenge is the quantification of theo¬ 
retical uncertainties [19]. Calculations in oxygen and cal¬ 
cium isotopes based on nuclear forces derived from chi¬ 
ral effective field theory (EFT) [20, 21] suggest that the 
uncertainties from the many-body methods are well con¬ 
trolled [4, 5, 11, 13, 22]. Therefore, uncertainties in the 
input Hamiltonian, such as truncations in the chiral EFT 
expansion or uncertainties in the low-energy couplings, 
likely remain the dominant source of uncertainty. Note 
that recently, first studies of the statistical uncertain¬ 
ties from numerically optimized chiral forces [23, 24] and 
to quantify correlations between chiral EFT couplings 
[25, 26] have been performed. 

In this work we investigate all sd-shell nuclei based 
on chiral two-nucleon (NN) and 3N interactions with 
realistic saturation properties. We derive microscopic 
valence-space Hamiltonians, which we diagonalize to ob¬ 
tain ground-state energies, two-neutron and two-proton 
separation energies, and first excited 2f energies. By 
varying the resolution scale in nuclear forces and the low- 
energy 3N couplings, we provide theoretical uncertainty 
estimates for these observables. In addition we also ex¬ 
plore the uncertainty associated due to the many-body 
calculations. We find that the resulting energy ranges, 


which are in good agreement with experimental data 
throughout the region, are mainly driven by uncertainties 
in the Hamiltonian. 

Nuclear interactions. At the NN level, we start 
from the next-to-next-to-next-to-leading order (N 3 LO) 
500 MeV potential of Entern and Machleidt (EM) [27]. 
We then use the similarity renormalization group 
(SRG) [28, 29] to evolve this interaction to a series of low- 
resolution scales A nn = 1.8,2.0, 2.2fm _1 . Taking chiral 
EFT as a general low-momentum basis, we combine each 
SRG-evolved NN interaction with the leading N 2 LO 3N 
forces [30, 31], where the Ci couplings in the two-pion- 
exchange 3N interaction are taken consistently with the 
NN interaction: C\ = — 0.81 GeV -1 , c 3 = —3.2 GeV -1 , 
C\ = 5.4 GeV -1 . We also vary independently the 3N 
cutoff A 3 n = 2.0,2.5fm -1 for the Ann = 2.0 frn -1 in¬ 
teraction. In addition, to probe uncertainties in the c, 
couplings, we use 3N forces with the Cj values obtained 
from the Nijmegen NN partial wave analysis (PWA): c\ = 
-0.76 GeV -1 , c 3 = -4.78 GeV -1 , c 4 = 3.96 GeV -1 [32] 
for the Ann/A 3 n = 2.0/2.0fm -1 interaction. For all 
Hamiltonians, the low-energy couplings cd, Ce in the 3N 
one-pion-exclrange and 3N contact interaction have been 
fit to the 3 H binding energy and 4 He charge radius us¬ 
ing Faddeev- and Fadeev-Yakubowsky calculations [15]. 
With this set of five NN+3N interactions, which predict 
nuclear saturation properties within uncertainties [15], 
we explore all sd-shell nuclei. The different interactions 
are denoted as 1.8/2.0 (EM), 2.0/2.0 (EM), 2.2/2.0 (EM), 
2.0/2.5 (EM), 2.0/2.0frn -1 (EM+PWA), where the la¬ 
beling indicates Ann/A 3 n and the c,; couplings used. 

Microscopic valence-space Hamiltonians. Based on 
the above interactions, we construct an effective valence- 
space Hamiltonian, H c g = + Peff, where £i de¬ 

note the single-particle energies (SPEs), and V/ff is the 
effective two-body interaction for valence nucleons. With 
H e ff we perform valence shell-model calculations, where 
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the many-body problem is solved exactly for the particles 
in the valence space on top of a closed core. 

Recently, nonperturbative methods for calculating H e s 
in medium-mass nuclei have been developed [33-36] and 
applied to the oxygen and very recently fluorine [37] iso¬ 
topes. In this work, we use many-body perturbation the¬ 
ory (MBPT) [38], which provides a diagrammatic order- 
by-order expansion for SPEs and V e s, to take into ac¬ 
count many-body processes outside the valence space. 
At third order, MBPT based on RG/SRG-evolved in¬ 
teractions with low cutoffs shows a reasonable order-by¬ 
order convergence for SPEs and V e g in medium-mass nu¬ 
clei [39]. To explore uncertainties associated with MBPT, 
we study valence-shell Hamiltonians obtained at second- 
and third-order MBPT. In contrast to phenomenological 
interactions, such as USDA/B [40] that fit both SPEs and 
V eS to experimental sd-shell data, our results are with¬ 
out adjustments. Therefore, we do not expect to reach 
accuracies comparable to the best fit USD interactions. 

Studies of oxygen [41, 42] and calcium [8, 10, 39] iso¬ 
topes showed that extending the valence space beyond 
one major shell provides additional binding and can lead 
to improvements for neutron-rich systems. However since 
our main interest is to perform a comprehensive study 
of sd-shell nuclei, estimating the theoretical uncertainty 
associated with the initial Hamiltonians, we limit our va¬ 
lence space to the sd shell (d 5 / 2 , d 3 / 2 , and s 3 / 2 proton and 
neutron single-particle orbitals on top of a 16 0 core). We 
work in a harmonic-oscillator basis with Juo = 13.53 MeV, 
appropriate for the sd shell, and scale all matrix elements 
of U e ff and bound SPEs by A -1 / 3 to correct for the in¬ 
crease in nuclear size. For all Ann considered, the cal¬ 
culations are converged in a basis consisting of 13 major 
shells for NN forces [39]. For 3N forces, we allow a to¬ 
tal energy of the three single-particle states up to 12 fux> 
(E 3 max = 12) in a basis of 13 major shells. 

Results. Figure 1 shows the ground-state energies 
of magnesium and chlorine isotopes compared to the 
Atomic Mass Evaluation (AME 2012) [43]. The second- 
and third-order MBPT results are represented by the 
blue, darker and the cyan, lighter bands, respectively, 
where the width of each band is spanned by the five dif¬ 
ferent NN+3N interactions considered. 

The experimental ground-state energies for magne¬ 
sium isotopes are generally within our uncertainty band, 
with neutron-rich isotopes at the lower side. Only the 
most neutron-rich isotopes are underbound in our calcu¬ 
lations. On the other hand, the ground-state energies of 
all chlorine isotopes are in good agreement with our un¬ 
certainty band, and they are typically within the lower 
side of the the third-order MBPT band, defined by the 
Ann/A 3 n = 1.8/2.0, 2.0/2.0, 2.2/2.0, 2.0/2.5 fm" 1 inter¬ 
actions. In general we find better agreement between our 
results and experiment for the isotopic chains of heavier 
elements, which suggests a somewhat too weak neutron- 
neutron interaction in our sd-shell calculations. This was 
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FIG. 1. Ground-state energies of the magnesium (top) and 
chlorine isotopes (bottom panel) relative to 16 O at second 
(blue, darker band) and third order (cyan, lighter band) in 
MBPT and compared to the Atomic Mass Evaluation (AME 
2012) [43]. The uncertainty bands are spanned by the five 
different NN+3N interactions (see text for details). The or¬ 
dering in the legend is with decreasing ground-state energies. 


also observed in Refs. [41, 42] for the oxygen isotopes. 

The estimated uncertainties in calculated ground-state 
energies are dominated by the different input Hamilto¬ 
nians. Specifically, the resolution-scale dependence by 
varying Ann from 1.8 — 2.2fm~ 1 with A 3 n = 2.0fnU 1 
is somewhat larger than the A 3 n dependence from 2.0 — 
2.5 fnU 1 for Ann = 2.0 fm -1 . This results in a combined 
resolution-scale dependence of approximately 1.0 MeV 
per valence particle in 32 Mg and 37 C1. When also in¬ 
cluding the PWA Ci values in 3N forces, the uncertainty 
roughly doubles to about 2.0 MeV per valence parti¬ 
cle. For both second- and third-order MBPT bands, the 
2.0/2.0fnU 1 (EM+PWA) interactions generally define 
the least bound calculations (for ground-state energies, 
the upper end of the bands). 

The difference between second- and third-order MBPT 
results is relatively small compared to the width of each 
band, indicating a reasonable, but still incomplete, con¬ 
vergence of the MBPT approach in this region. For mag¬ 
nesium and chlorine, third-order results are more bound 
because of more attractive proton-neutron interactions, 
whereas for oxygen (not shown), second-order results are 
more bound than at third order mainly due to the neu¬ 
tron single-particle energies. When the uncertainty as¬ 
sociated to the MBPT is also included, the total uncer¬ 
tainty increases to 2.1 MeV and 2.8 MeV per valence par¬ 
ticle in 32 Mg and 37 C1, respectively. 
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FIG. 2. Uncertainty estimates for the two-neutron separation energies S^n of sd-shell isotopic chains at second (blue, darker 
band) and third order (cyan, lighter band) in MBPT and compared to the Atomic Mass Evaluation (AME 2012) [43]. 


Figure 2 compares theoretical and experimental two- 
neutron separation energies S 2 n for all isotopic chains 
from oxygen to calcium (Z = 8 — 20). The theoretical 
calculations describe the overall experimental trends rea¬ 
sonably well, but in general our uncertainty bands under¬ 
estimate the empirical values. This is especially the case 
in lighter elements and for the most neutron-rich nuclei 
for all isotopic chains. This is probably related to the 
underbinding of the sd-shell calculations when valence 
neutron-neutron interactions are dominant. We also note 
that around N = 20, the ground states of 29,30 Ne [44, 45], 
30,31 Na [46, 47], and 31,32 Mg [48, 49] are dominated by 
deformed configurations not captured in our sd-shell cal¬ 
culations (this is the so-called island of inversion). Conse¬ 
quently, our bands do not reproduce the change in slope 
of S' 2 n around N = 20 for Ne, Na, or Mg. 

Similar to the ground-state energies, the dominant un¬ 
certainties arise from the different Hamiltonians, with 
smaller differences between second- and third-order 
MBPT results. Typical the uncertainty range for S 2 n is ~ 
5 MeV. The exception are N < Z isotopes, more visible 
in heavier elements, where the difference between second- 
and third-order results is comparable to the uncertainty 
between input Hamiltonians, due to too weak proton- 
neutron interactions at second-order MBPT, adding up 
to a total uncertainty of ~ 10 MeV. 

In Fig. 3 we show the two-proton separation energy S^p 
for all isotonic chains from TV = 8 to N = 20. Our re¬ 


sults agree very well with experiment in all cases, and re¬ 
markably most experimental values fall within the third- 
order MBPT band. Only in few proton-deficient and 
very proton-rich nuclei do experimental lie within the 
second-order band. Since there are fewer proton-rich nu¬ 
clei known experimentally than neutron-rich nuclei, S^p 
are in general informative about proton-neutron interac¬ 
tions. The much better agreement in S 2 P than for S 2 n 
compared to experiment suggest that the different Hamil¬ 
tonians considered capture better (mostly isoscalar) 
proton-neutron interactions than neutron-neutron inter¬ 
actions. Again, the sensitivity to the input Hamiltonians 
dominates the theoretical S 2 P uncertainties (with a simi¬ 
lar range of ~ 5 MeV), except for proton-deficient nuclei 
where the MBPT uncertainty is comparable (with a total 
uncertainty of ~ 10MeV). 

Finally in Fig. 4, the calculated first excited 2+ en¬ 
ergies are compared to experimental data for all even- 
even sd-shell isotopes. The spread of the uncertainty 
band is typically smaller than ~ 500 keV, with generally 
reasonable agreement to experiment. However, in the 
cases with high-lying 2]*~ states, indicative of shell clo¬ 
sures ( 22 0, 24 0, 22 Si, 34 Si, 34 Ca), the uncertainty can be 
as large as ~ 1 MeV. This means that, while our bands in 
general predict shell closures consistently, the actual exci¬ 
tation of the state is very sensitive to the details of the 
input Hamiltonian. The width of the uncertainty band is 
mostly due to the 2.0/2.0fm _1 (EM+PWA) interaction, 
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FIG. 3. Uncertainty estimates for two-proton separation energies S 2p of sd-shell isotonic chains at second (blue, darker band) 
and third order (cyan, lighter band) in MBPT and compared to the Atomic Mass Evaluation (AME 2012) [43]. 


which is also responsible for the unusually large uncer¬ 
tainty band in 36 Ar. In general, the second- and third- 
order MBPT bands mostly overlap, except for TV ~ Z ar¬ 
gon and calcium isotopes, where only third-order MBPT 
results are in agreement to experiment. Similar to the 
case of S 2n , we also note that the 2j states within the 
island of inversion, the TV = 20 isotopes 30 Ne [45] and 
32 Mg [49], are deformed, and their relatively low exci¬ 
tation energies cannot be well described in our sd-shell 
calculations. 

Summary. We have presented a comprehensive study 
of ground-state energies, S 2n , S 2p , and first excited 2± en¬ 
ergies for all sd-shell nuclei: isotopic chains from oxygen 
to calcium and isotonic chains from TV = 8 to TV = 20. 
This is based on NN+3N Hamiltonians that have been 
fitted only to A = 3,4 nuclei that predict realistic sat¬ 
uration properties of nuclear matter, without additional 
adjustments. We have focused on estimating the theoret¬ 
ical uncertainties due to the different input Hamiltonians 
and associated with the many-body calculations. We find 
reasonable agreement to experimental data, especially in 
nuclei dominated by valence proton-neutron interactions. 
For neutron-rich systems, calculations in extended va¬ 
lence spaces are needed, due to too weak neutron-neutron 
interactions in the sd, shell. 

Generally, we find that the estimated theoretical un¬ 
certainties are dominated by differences in the NN+3N 
interactions. While the present first study is limited to 


the EM 500 MeV N 3 LO NN potential, it will be im¬ 
portant to perform more comprehensive studies. These 
efforts will incorporate the exploration of different fitting 
procedures for the 3N low-energy constants cd,Ce, dif¬ 
ferent regulator forms and cutoff values for NN and 3N 
interactions as well as order-by-order convergence stud¬ 
ies in the chiral EFT expansion [51-54]. For improved 
studies of the resolution scale dependence, we plan to 
perform calculations based on consistently SRG-evolved 
NN and 3N interactions [55, 56]. Improving our uncer¬ 
tainty estimates due the many-body calculation is more 
challenging, because corrections to MBPT (i.e., fourth- 
order contributions) are not attainable at present. How¬ 
ever, nonperturbative methods for valence-space Hamil¬ 
tonians can provide a controlled framework to assess the 
many-body approximations [33-36]. 
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FIG. 4. Uncertainty estimates for excitation energies of the first 2f states in even-even sd-shell isotopes at second (blue, darker 
band) and third order (cyan, lighter band) and compared to experimental data from ENSDF [50]. 
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